Local entanglement generation in the adiabatic regime 
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We study entanglement generation in a pair of qubits interacting with an initially correlated 
system. Using time independent perturbation theory and the adiabatic theorem, we show conditions 
under which the qubits become entangled as the joint system evolves into the ground state of the 
interacting theory. We then apply these results to the case of qubits interacting with a scalar 
quantum field. We study three different variations of this setup; a quantum field subject to Dirichlet 
boundary conditions, a quantum field interacting with a classical potential and a quantum field that 
starts in a thermal state. 
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PACS numbers: 03.67.Bg, 03.70.+k 

I. INTRODUCTION 

Quantum entanglement is widely believed to be the 
distinguishing resource of quantum computers. For this 
reason it is crucial that we fully understand how entan- 
glement can be generated and how it evolves. Various 
studies already show that the time evolution of entan- 
glement in open systems is often highly non trivial, see 
e.g. [l|. This complicated time evolution often prevents 
us from studying general features of entanglement dy- 
namics and forces us to focus on particular examples. In 
this work, we focus our attention on one particular time 
evolution scenario, namely, the adiabatic evolution of the 
ground state. This allows us to show that adiabatic evo- 
lution can naturally generate entanglement in a pair of 
qubits interacting with a correlated system. 

One of the main applications of this setup is the ex- 
traction of entanglement from the vacuum. Indeed, it was 
shown in |2( that qubits interacting adiabatically with a 
rclativistic quantum field in the vacuum state can get en- 
tangled in a renewable fashion. This result is one of the 
many new features that arise as a consequence of special 
relativity considerations in quantum information theory, 
see e.g. 0-0]. We follow-up on this work by studying 
modifications of this setup and analyzing whether they 
enhance or degrade the entanglement generated in the 
qubits. 

The paper is organized as follows. In Sec. II we present 
the general framework and calculate explicitly the entan- 
glement contained in a pair of qubits in the ground state 
of a weakly interacting theory. We also discuss how this 
entanglement can be generated with an adiabatic switch 
on of the interaction Hamiltonian. In Sec. Ill we use 
the tools previously developed to show that the entangle- 
ment available in a quantum field theory with Dirichlet 
boundary conditions is degraded. In Sec. IV we consider 
a quantum field weakly interacting with a classical field 
and show that depending on the type of interaction it 



can either enhance or degrade the entanglement gener- 
ated in the pair of qubits. In Sec. V we study the case 
of a quantum field that starts in a thermal state instead 
of the vacuum and show how the entanglement decreases 
as the temperature of the system increases. 

We work in the natural units K = c = 1. Wherever 
necessary to avoid ambiguity we will denote operators O 
or states |t/»), corresponding to the Hilbcrt space Hy\ by 
a superscript (j), for example, O^' and |'0 )■ Orders in 
perturbation theory will be denoted by a subscript (j), so 
for example we could have P = P( ) + Pm + 0(a 2 ). We 
conveniently work in the Schrodinger picture of quantum 
mechanics. 



II. GENERAL FRAMEWORK 
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FIG. 1: Qubits A\ and A2 interacting with system B. 

The system we study is illustrated in Fig. ((T|) which 
describes two localized qubits, A\ and A2, unitarily and 
locally interacting with a correlated system, B. Let the 
Hamiltonian of the free theory be of the form Hq = Hb + 
Ha- More precisely, let 

2 

H = ^Eklk^) (k {B) \+J2[^ E s + AE ) 
k j=i 

x\e^)(e^\ + E g \g^)(g^\]. (1) 

The assumption that both two-level systems are identi- 
cal allows setting E g = 0. In addition, for system B, 
we set E Q — and denote the corresponding eigenstate 
by Ig^)- Thus, the ground state of the free theory is 



\g {B \g {A) ) (where \g^) := \g {Al) , g iM) )) and we have 
Hq \g^ B ',g^ A ') = 0. It is convenient to choose the basis 
{|s (j4) )} (for n (Al) <g> 'H (A2) ) coinciding with the eigen- 
states of Ha, that is 

{| S ^>} = {|l( A )) = |e^,e^),|2^) = |e^ jfl (>i a )) j 
|3< A >) = |^ l) J e^)),|4^)) = |^), 5 ^))}. (2) 

Then we have H \k^ B \ s^) = (E k + e s ) \k^ B \ s^} 
where t\ = 2AE, e 2 = 63 = AE and 64 = 0. 

Let the local interaction between the qubits A\ and A 2 
and system B be described by a Hamiltonian of the form 



Hi, 



Hba x + Hba 2 - 



(3) 



For simplicity we assume that 

(g {B \g {A) \H m t\g iB \g {A) ) = 0. Furthermore, we 
assume that the interaction described by H int is weak, 
so we can resort to perturbation theory to find the 
ground state of the interacting theory |f2), that is, 
(Hq + H int ) |fi) = Eq |fi). Thus, up to second order, we 
may write |fi) = |fi(o)) + |^(i)) + |^(2)) + 0(a 3 ) where 
a is a small parameter that sets the scale of Hi nt . We 
have |f2( )) = \g^ B \g^) and using time-independent 
perturbation theory [8[ one can easily find explicit 
expressions for |f2m) and l^m)- The density matrix 
describing the joint system AB after the interaction 
reads p := \Q) (Cl\ = /9 (0 ) + /3(i) + p{2) + 0(a 3 ) where 
P(o) = l^(o)) (fyo)|, P(i) = (|tyo))(tyi)| + h.c.) 
and p(2) = |0 (1) ) (0 (1) | + (|fi( )) (fi( 2 )| + /i.e.)- Using 
these expressions, we readily determine the reduced 
density matrix describing system A\A 2 when the 
joint system AB is in the ground state of the inter- 
acting theory. Thus, up to second order, we have 
p A ■= Ti B (p) = PA(o) + PA(i) + PA(2) + 0(a 3 ) with 

PA { v) = \g {A) ){g (A) l 



1/ (B) 8 {A)\ H . I (B) JA)\ 

PAH) = -E - } \s iA) )(9 (A) 

+h.c, 

^(2) - -|5 l X.9 J IE (^ +es)2 



(4) 
2 



-E 

k,s,r 



(g( B \g( A )\H mt \k( B \r( A )) 

(E k + Cr) 



/ fc (B) (A) I u-. I (B) Q (A)\ 

V^ if I -"»"« Iff iff / 1 (A)\ / (A) I 

(E'fe +e s ) 



E 



(ffW.^l^tlfc^.r^) 



(£ fc + e r ) |S K9 ' 



+h.c. 



(5) 



Here the sums J^ fc s run over all the values of {k, s} 
except those for which any denominator vanishes. For 
simplicity let us now focus on an important class of 
local interactions, namely HsA k = afc(|e' Afc ^) (<r fc '| + 
| e (A fc )^ 5 (A fc )Q F (B) ; for fc = (1,2). The nonvanishing 
matrix elements are then given by 



Pi 



Pi ■ = 



>?E 

k 

*E 



E 



F 



ai«2 E 

k 



(g^\¥ {B) \k^)\ 2 

(E k + aey- 

(g^\¥ {B) \k^)\ 2 
(E k + AE) 2 

(g( B )\¥ (B) \k^)(k^\¥ (B) \g^) 



(6a) 
(6b) 



E 



(E k + AE) 2 
(g^\¥[ B} \k^)(k^\¥ { 2 B} \g^) 



(6c) 



AE(E k + AE) 



(6d) 



PA 



'0 F* 

Pi E* 

E P 2 

y F 1-Pi 



+ 0(a 4 ) (7) 



P7 



We note that the above matrix elements may be 



written as 



7WI 



7( b )im1 b ) 



fy*)) 



A 19 

1 {B) and Nf 
,a k ¥[ B) and Nf 



2AB 

tors M. k "' and N)."' for fc = (1,2) are defined as ivu fc 



,a k ¥^ 



P-2 = 

g {BS >) and 

The opera- 
te _ 

k ~ 

From these 



Hb+AB"* 1 k allu iN fc _ ^/H B +AE 

expressions one easily proves that P1P2 > \E\ 2 , thus 
guaranteeing the positivity of pa ( up to second order 
in a). In order to quantify the degree of entanglement 
in system A\A 2 we make use of the negativity N(pa), 
defined as twice the absolute value of the negative eigen- 
value of p TA i Q. In our particular case it reads 



N{ PA ) = max (ytPi-p^+^l 2 - P x - P 2 , 
+0(a 4 ). 



(8) 



To generate this entanglement in the pair of qubits, we 
need to prepare the state of system AB in the ground 
state of the interacting theory |0). To do this, we assume 
that the state of system AB can easily be prepared in the 
ground state of the free theory \g^g^ A ')- Moreover, we 
assume that the interaction Hamiltonian can be switched 
on with a switching function r)(t) such that H(t) = Hq + 
T}(t)H int where rj(t < U) = and r](t > t. h + At) = 1. 
If the interaction between the qubits and B is switched 
on adiabatically, then the evolution of the joint system 
is given by \g B ,g ) — > |fi(i)) where |f2(i)) is the ground 
state of H(t) . According to the validity condition for 



adiabatic behavior [10, llll , we need at first order 



,-/m ■ (E k + AEV 
max \mt)\ <C min 7=7 

^ Ul<s (fl) |F ^^L 



(9) 



Therefore, if this condition holds for some choice of rj(t), 
then the ground state |0) of H(ti + At) can easily be 
reached in a time scale of At ~ 1/ maxt |?y|. 



A. Example: Qubits interacting with a scalar 
quantum field 

As an application of this formalism, we consider qubits 
interacting locally with a smeared portion of a quantum 
scalar field <j>(r) of mass m. This effectively models an 
atom interacting with a quantum field like the quantum 
electromagnetic field. This example was first explicitly 
considered in [2|. The operators Fjj. are then 



F 



(23) 



d 3 rf k (r)0(r). 



(10) 



and for simplicity we choose a± = a 2 = a and /2(p) = 
/i(r — d) such that the distance between the two qubits 
is d. Note that the smearing functions /& (r) describe the 
effective size AX of the qubits. In the limit -^ — > 00 the 
introduction of the smearing functions /& (r) is equivalent 
to the introduction of a cut-off A ~ -^ in momentum 
space. Therefore, we shall always replace these smearing 
functions with a momentum cut-off and set 



r(B) 



(Tit 



From Eq. (foaj) . (|6d|) and © we recover the results of 



a 2 f 1/AX 



Q 

4^ 



£ p (£ p + A£) 2 



4 / AX psin(pd) 

P £p(i5~fAB)(ABdy 



AA = 2max(|P|-P,0)+O(a 4 ) 



(11) 



(12) 
(13) 



where E p = \fp^ 
in the limits dAE 



easily 



show 



2W 1 LLLa ^ \ 2dAE 

m we have M 



that if 
-In 



AE 

1 > 



AEAX 

max 



Using these equations 

and dm — > one can 

3> m we have M ~ 

, 0) and similarly if AE <C 

ln(^r),0). More- 



2(/A£ 



Z7T- \ AUjLXTj \ rilL\s\ I ' / r— I 

over, using Eq. ([9]) with Eq. ([T0|) it was show in |2| that 
adiabatic evolution is possible in principle. In fact, in 
order to have a very small error in the ground state neg- 
ativity at the end of the time evolution we roughly need 
max t \f\[t)\ -C AE. Thus, if we follow that prescription 
we can adiabatically switch on the interaction and keep 
all a 2 contributions in Eq. (fT3"|) intact. 



III. QUANTUM FIELD WITH BOUNDARY 
CONDITIONS 



In this section we follow-up on the previous example 
by considering the case where the scalar quantum field is 
subject to Dirichlet boundary conditions. Such scenarios 
naturally arise when describing electromagnetic waves in- 
teracting with perfect conductors and have been exten- 
sively studied in the context of the Casimir effect [13 |. 
Here, our goal is to investigate whether the presence of 
boundary conditions augments or reduces the amount of 
entanglement generated in system A\A 2 - For simplicity 
we only consider a massless field. In this case, the field 
operator is expanded in terms of creation and annihila- 
tion operators as 

4>{r) = J2 -i (wW + 4*4^)) ( 14 ) 



Op , at, 



Up.p' 



(15) 



where the up{r) are solutions of Helmholtz equation (A + 
\p\ 2 )up{r) = satisfying the boundary conditions. The 
matrix elements Pi, Pi and F arc easily expressed in 
terms of the mode functions up(f). From Eq. (|6a|) . (J6b| 
and (I6dl) wc obtain: 



Pi 



P> 



*?£^ 



v(^i)l 5 



- 2|pl (|pl + AEf 



a ^E^ 



M^)| s 



- 2|pl (|pl + AE) 2 



F = aia^SR 



y- 



(Mj?(fi)ut(r 2 )) 



^2\p\(AE)(\p\+AE) 



(16a) 
(16b) 

(16c) 



Let us consider the scenario in which the field 4>(r) satis- 
fies the Dirichlet boundary conditions tfi(±kf,y, z) = 0. 
In addition, we temporarily impose periodic boundary 
conditions on the y-z plane, i.e. 4>{x, y + L y , z + L z ) = 
<j>(x,y,z). Under these assumptions, the mode functions 
Uf(r) read 



Px [x + — 



y L y L z 



(17) 



and 



where p\\ = (0, 2*2*, 2%±) and p x = ^. Here 
n z assume the values 0, ±1, ±2, . . . whereas n x = 1,2, . . .. 
Note that in Eq. (|16a[) , (|16b>|) and (|16c[) we need to deter- 
mine sums of the form ^_. c(jT)up(fi)u*Jf2). In the limit 
00, these sums take the form 



L y —> 00 and L 



p 



c{p)Up{ri)u}{r 2 ) 



E 



d 2 P\\ c(n,p||) e ^ r(ft _ ft) 



(2tt) 2 2L 

-(xi-x 2 ) _ 



~-(x 1 +x 2 +L 3 



(18) 



Making use of Poisson summation formula 



above sum as 



^c(p)u p -(fi)u;-.(r 2 ) = Y^ 

p nEZ 



= Enez 3( x ~ n ) one can rewr ite the 



d 3 p 

(2^)3 



>m 



,ipR n _ e ipR n 



(19) 



where R n = (xi - x 2 + 2nL x ,yi - 2/2,21 - z 2 ) and 
i?„ = (xi + x 2 + (2n + l)i a! yi - 2/2,21 - z 2 ). From 
the above equations one can determine the entanglement 
in A\A 2 for arbitrary positions of the qubits. We will 
however limit our discussion to two symmetric configu- 
rations. Let us first consider a symmetric configuration 
such that the qubits are located at ffc = (±|,0, 0) with 
[d < L x ). Assuming a\ = a 2 — a and making use of Eq. 
(jT6aJ), ()16cj) and (fT9|) we arrive at the following expres- 
sions: 



P 



F 



|p|<l/AX 



x ( e lp * 2n 

- 2 E 

nGZ 

(< 



d 3 p 



(27r)3 2|p|(|p| + A£;)2 



ip^(d+(2n+l)L a; ) 



d 3 p 



x e 



|p|<l/AX ( 27r ) 3 2 
,ip x (d+2nL !c ) _ ip a ,(2n+l)i a 



AE)AE 



(20) 



(21) 



Note that the free space situation (i.e. in the absence 
of boundary conditions) may be recovered by taking the 
limit L x — > 00. Indeed, in the regime L x 3> d, Eq. (|2"0)l 
and dHJ reduce to Eq. (JTTJ and ([11]) (with m = 0). It is 
convenient to express Eq. (|2T)|) and ([2"Tjl in terms of the 
dimensionless quantities \q\ = \p\L x , e = dAE, 7 = 77- 

and A = d/ AX. After simple manipulations we obtain 



A/7 



1 



47r 2 /_ d( i (, 



E 



(<7 + e/7) : 
sin(2n(7) sin ((2n + 7 + 1) q) 



F 



n 

a 7 



2» 

A/7 

dq- 



2n + 7 + 1 



(22) 



1 



4"7T 2 £ ,/ "* (7 + e/7 

sin((2n + 7) </) sin ((2n + 1) q) 



E 



2?i + 7 



2?7 + 1 



(23) 



Finally, by means of the formula [li 



Esin ((2n + a) a) tt /. . 7ra\ , „ 

2n + a - 2^n(f) ^ (^ + J > t) ^ 



for m7r < q < (m + l)ir we reduce Eq. (J22J) and Eq. ([23 
to the simpler form 



P 



a 

8^- 



M„ 



8^2 Z^ 



2m + 1 



07 
8rc 



a/ 



' (m + — )(m+ — + 1 

sin ((2m + 1) (7 + 1) f ) 

(2m+l)sin(( 7 +l)f) 

l\ 



(25) 



E^ 



77T 



m=0 

sin((2m + l)^ 



771 + — 

77T 



Sill 



' , 77T'\ 



(-1)' 



(26) 



where M max w A/(ttj). Note that in the limit 7 — >• 
1 the above expressions vanish in accordance with the 
boundary conditions. Consequently, the entanglement in 
the qubits should vanish as L x — > d. Numerical results 
for the entanglement generated in system A\A 2 versus 
7 = 77- are presented in Fig. @. Another particularly 




FIG. 2: K = ^rN as a function of 7 = 7^ £ (0.01, 1) with 
A = d/AX — 10 3 . The upper curve (red) corresponds to 
e — 0.015, the middle curve (blue) to £ = 0.02 and lower 
curve to e = 0.03. 

interesting case is that where the qubits are located at 
ffe = (0,±|,0). Making use of equations (f!9| we obtain 
expressions analogous to (|22[) and ([2U)l . They read 



P = 



4^E / ^ 



(g + e/7) 2 



F = 



a 7 

47r 2 e 



sin(2ng) sin((2n + l)g) 

277 

E 



«e: 



A/7 



dg 



2?7+l 

1 

? + e/7 



(27) 



3in( A /(2n) 2 + 7 2 g) sin(y / (2n + l) 2 +7 2 g) 



V(2nF+l 2 



V(2n+l) 2 + 7 2 



(28) 




FIG. 3: K = l^rN as a function of 7 = -±- £ (0.01, 2) with 
A = d/AX = 10 3 . The upper curve (red) corresponds to 
e = 0.015, the middle curve (blue) to e = 0.02 and lower 
curve toe = 0.03. 



Clearly, in this case the boundary conditions do not imply 
that matrix elements P and F should vanish as L x — > d. 
Numerical results for this configuration are presented in 
Fig. ([3]). Thus, both graphs indicate that the entangle- 
ment generated in the pair of qubits reduces monoton- 
ically as the separation L x decreases. Note that in the 
regime L x ^> d, the orientation of the qubits relative to 
the planes x = ±L x /2 becomes irrelevant and as a con- 
sequence the negativity values coincide for the two cases 
considered. 



IV. QUANTUM FIELD INTERACTING WITH 
A CLASSICAL POTENTIAL 



In this section we study entanglement generation in 
the qubits when system B is either self-interacting or 



J 



interacting with an external classical system. To do so, 
we consider the Hamiltonian 



H = Hq + H^ 



(29) 



where H mt = H int 

^2 _ nJA^\ iJA k ) 



XV ^ and as usual Hi nt = 

ELi M\e (Ak) ) (g (Ak) \ + \e {Ak) ) {g {Ak) \W { k ] '. Here V^ 
is a potential acting solely on system B. Throughout this 
section, we assume that 



(g B \V^\g B )=0. 



(30) 



Clearly, the presence of the potential V^ modifies the 
density matrix pa ■ Let us denote the new reduced density 
matrix by pa '■= Pa + &P\ where 8p\ contains all the 
contributions coming from the potential V^ B \ Following 
similar steps to those in Sec. [Til we a Pply perturbation 
theory to find the second order term (containing terms 
of the form CYfcA) 



5p\(2) = A 



E 

k,s 



I 1 



Ek^s 



(k^\ s ^\HM B \g {A) ) 



x { g ^\V^\k^) + {g( B \ S ^\H mt \k^\s^) 



x (k^\V^\g^))\ S ^)( 9 ^\ + h.c. 



(31) 



We note that for potentials built out of even powers of 
the field </>(r), the above expression vanishes. In order 
to include this class of potentials into our framework, we 
need to include third order corrections. Making use of 
the condition Eq. (|30]l . we obtain after some algebraic 
manipulations the third order correction to the reduced 
density matrix. It reads: 



PA(3) 



s,k,r,j,l 



(gC^.aWlfl-^lfc^.r^) <fc^,rW|g iTtt |j(^),Z^)> 0-t B ),ZW|g iTlt | g ( B ), g ( A )) (A) (A) 



(E k + e s )(E k +er)(E j + ei) ^ M * l + 



E' 

s,k,r,j 



(.g^),^)!^!^),,^))^^),^)!^!^)^-^))^)^-^)!^^),^^) 



(E k + e s )(E r + ej) 



1 



1 



Ek +e s E r + £j 



M)\ 



(g {A) \ 



r 



(32) 



Here, we note that the matrix pa keeps its original form relevant modifications of the matrix elements are given 
(as in (O). In other words, the corrections do not gen- 
erate new non- vanishing entries in the matrix ([?])• The 



by 



theorem Il5l we obtain 



SPi> 



-\alJ2 



(gW^^kW) (kW\VW\sW 



(E k + AE) 
1 1 



(E S +AE) \E k + AE ' E S +AE 

(E k + AEf 
( s (B)|^(B)| g (B) ) ( g (B)|y(B)| s (5)) 

E S Eg 

( S W\w[ B) \kW)(kW\F[ B) \gW) 



5F\ = — \a.ia2$t 2 



{E k + AEf 

(g^\¥[ B) \k^) 



(33) 



-s . A: 



AE(E k + AE) 



(kW\VW\sW) (sW\VJ B) \gW) 
(E s + AE) 

(g(B)\ V (B)\ k (B) ){k (B)\ ¥ (B)\ siB)) 

AEE k 



(E s + AE) AE 

(kW\w[ B) \sW)( S W\F { 2 B) \gW) 
E k {E s + AE) 



(34) 



Naturally, 8P2X may be obtained from SP±\ upon replac- 
ing F[ B) by ¥ {B) and a x by a 2 in Eq. ([33) . Note that 
when y( B ) = iJ B then 5P\\ and <5F^, obtained from the 
above expressions, coincide with the first order term ap- 
pearing in the Taylor expansion of Eq. (pa)) and (|6d[) . 
That is, we have Pi\E k ^(i+\)E k ~> -Pi|.E fc + <5-Pla plus an 
analogous relation for F. 

Let us now follow up on the proposal by Achim Kcmpf 
14 1 to study the case of qubits interacting with a quan- 
tum scalar field which is interacting with a classical po- 
tential. We model this by choosing 



V {B) = / d 3 rV (F) : 2 (f) 



(35) 



The normal ordering :: [15| automatically guarantees that 
condition (|30[) is satisfied and it renders the matrix ele- 
ments of Eq. ([35)) and ([34]) finite. This model can be seen 
as an analog of QED where the electromagnetic field is 
in a coherent state and therefore can be treated classi- 
cally. For this reason the model is very similar to po- 
tential problems in non-relativistic quantum mechanics. 
We now proceed to compute the corrections 6Pi\ and 
SFx given by Eq. (J33J and ([54"]). Making use of Wick's 



SP 1} 



Act 2 



d 3 Pi 



d 3 p 2 



2 4-iKi/Ax (^) 3/2 J\p 2 \<i/AX (2tt)3/ 2 

1 



V(P2-pi) -jfa-avft 

j -'pi '-'pi 

1 



£jp 1 + -&P2 



1 



(%+Afi) 2 (i% 2 +AF) 2 
1 
{E Pl + AE){Ep 2 + AE) 

1 1 



% 2 + AE Ef; 2 + AE 



SFx 



Xa 1 a 2 f d 3 pi 

~ 2A-B J Wi1<1/ ax (2^) 3 / 2 7| p - 2 |<i/a X (2tt) 

x V(P2 ~ Pi) c -i(g,-f a -pl.n) 



(36) 
d 3 p 2 



W2 



J^pi -^p 2 



1 



1/1 1 



% + £ft V % + A ^ %2 + AF 



(37) 



where V(p) := J ^ 7r w e I P' r V(r'). Here note that if 
we set V(r) = ^- then we are simply dealing with 

a Klein- Gordon Hamiltonian with mass •v/(l + A)m. 
The reader may check that in fact the above expres- 
sions reproduce the correct Taylor series expansions of 
(ED) and (J6HJ). That is, Pi|„^,. /ttttt,^ ~> P 
<5Pia and F 



m->-(- v /(l+A)m) 

-»■ PI™ + <5F A - We will 



| ™-kV( i +*)™) 

make use of this simple observation in the next subscc 
tion. 



A. Example: Spherically symmetric Gaussian 
potential. 

We now apply the above results to the situation where 
two identical detectors ai = a2 = a, located at r k = 
(±4,0, 0), are interacting with the massive scalar field 
(f>(r) coupled to the spherically symmetric Gaussian po- 
tential 



V[f) = V e 



(38) 



Its Fourier transform is easily found to be V(p) = 
j2 V ?3/2 (7 B e ~ lC ' B ^ ■ The axial symmetry of the problem 
may be exploited by making use of plane wave expansion 
into spherical harmonics. Thus, after some algebraic ma- 
nipulations, we obtain the useful identity 
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FIG. 4: K = ^-TV as a function of ma B - Here AE/m = 0.1, 
mAX = 10 -3 and md = 0.5. The middle dashed curve cor- 
responds to the case A = 0, the upper thick curve (red) cor- 
responds to Wo/m 2 — — Y^Q and the lower thin curve (blue) 
to Wo/m 2 — -j^q • Upper and lower dashed lines (orange) 
correspond to the values Kf(m — ¥ -^/l =F l/50m) with A = 0. 
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FIG. 5: K = ^-N as a function of ma B - Here AE/m = 
0.1, mAX = 10~ 3 and md = 0.9145. The thick curve (red) 
corresponds to Wo/m 2 — —j^ whereas the thin straight 
line (blue) corresponds to Wo/m 2 — -A^. The dashed line 
(orange) represents the value Af(m — >• a/1 — l/50m) with A = 
0. 
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(39) 



where J v (x) and I v (x) are Bessel functions [13]. The 
above expression facilitates considerably the numerical 
evaluation of the integrals in Eq. ([31)]) and ([3T]) . The 
effect of V(r) on the entanglement of system A1A2 is 
shown in Fig. (J4]) and ([5]). In Fig. (j4}, we chose a set 
of parameters AE/m, AEd and mAX such that the de- 
tectors are entangled as a result of their local interaction 
with the quantum field. We show the negativity of the 
qubits as a function of the width of the potential for the 
cases Vq > (repulsive potential) and Vq < (attrac- 
tive potential). In agreement with intuition, we observe 
that entanglement increases for the attractive potential 
whereas it decreases for the repulsive potential. More- 
over, note that when the width of the potential is much 
greater than the separation between the qubits (as ^ d) , 
they effectively experience a locally constant potential. 
As previously discussed, this situation is equivalent to a 



mass shift given by m 



l eff 



1 + 2A 



v 



In 



Fig. (01 we show the asymptotic values M(m — > m e ff) 
determined by Eq. ([II]), tfE]) and ([13]). Finally, in Fig. 
©; we chose a set of parameters AE/m, AEd and mAX 
such that Af=0 and |F| « P. In this case we see that the 
small correction to F and P coming from the attractive 
potential (Vo < 0) may induce entanglement in system 
AiA 2 when the potential is sufficiently wide. The results 
presented in this subsection simply reflect the fact that 
particle exchange between the detectors tends to be fa- 
vored by a central attractive potential and hindered by 
a repulsive potential. 



V. QUANTUM FIELD IN A THERMAL STATE 

In this section we consider the entanglement generated 
in the qubits when they are interacting with a scalar 
quantum field of mass m > AE known to be initially 
in a thermal state. In general, the adiabatic theorem 
and time- independent perturbation theory are tricky for 
thermal states because the energy levels are degenerate. 
However, in our case the interaction Hamiltonian does 
not couple states of equal energy when m > AE and 
therefore degeneracy poses no additional complications. 
In fact, note that even though system B is not initially 
in the ground state, one can easily verify using Eq. ([9]) 
that adiabatic evolution is still possible provided that 
m > AE. Let us first consider a general multi-particles 
state 0: 
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Replacing \g^) with |V> (S) ) in Eq. flBa} and J6d]) we find 
when ai = a2 = a'- 



P := P k = a 2 
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(41) 



F = q 2 K 



(AE - EpY t 

^ (V^| 0(rx) |^ + y) (^ + l p -| 0(r 2 ) IV;) 
(Ej!+AE)AE 



E 



(AE - E P )AE 



(42) 



Note that there is an additional term when we consider 
excited states, this additional term accounts for the pos- 
sibility that the held destroys an existing particle. Simple 
calculations show that: 



(V> + lp| 0(r) |V) = V 1 + n r 
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(44) 



Using these equations, one obtains in the continuum limit 
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(45) 



(46) 



Note that the vacuum situation may be recovered by tak- 
ing riff = 0. Indeed, when np = Eq. (|4"5)) and (|4"6")l 
reduce to Eq. (jTTJ) and (fT2|) . We now assume that the 
field starts in a thermal state of temperature T, that is 



(47) 



S /J% - i 



where /3 = -^. Using Eq. (J45]) and flU with Eq. (l47l) . 
we can numerically evaluate the negativity as a function 
of the temperature, see Fig. §§§. We can also find an 
explicit expression for the first order correction to the 
negativity in the low temperature regime ((3m 3> 1 and 
— tU 3> 1) under the assumption that m 3> AE. First 
note that in the limit AE/rn — >• 0, F is unchanged by 
the temperature. On the other hand, in the same limit 
P has a correction which we may denote as P(i) such that 
P = P( ) + P(i) where P( ) is given by Eq. (fTT|) and P^) 
reads 



P (D 



l/AX 



dp 2p 2 n p 
47T 2 (p 2 + m 2 ) 3 / 2 ' 



(48) 



In the low temperature regime we have np ss e -0 E p \ We 
thus see that np provides a smaller effective cut-off to the 
momentum integral than the momentum cut-off caused 
by the size of the qubits. Therefore, we can make the 
crude approximation 
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(49) 

As ex- 



such that we roughly have A/"~jv( — u - 2 .^ . 
pected, the negativity decreases with an increase of the 



temperature. This also implies the existence of a critical 
temperature, that is, it is possible to extract entangle- 
ment from the quantum field provided that its tempera- 
ture is below T c . This critical temperature is given by 



ksTc 



2m, 



(50) 



W 2 T 

\(2dSB- ln (^ix)) 

where W(rr) is the Lambert function [17| . 








FIG. 6: K = 2VA/" as a function of 6 = k B T/m with 
AE/m = 0.1 and mAX = 10~ 3 . The upper curve (red) 
correspond to e = AEd — 0.07, the middle curve (blue) to 
s = AEd = 0.075 and the lower curve to e = AEd = 0.08. 



VI. CONCLUSIONS 

In this paper we showed how entanglement dynamics 
can be studied in the adiabatic regime. In this regime, 
the time evolution of entanglement is relatively simple 
which allows us to study more exotic setups and still 
gather valuable insights on the behavior of entanglement. 
As an example we studied qubits interacting with a quan- 
tum field in non-trivial contexts and we arrived at the 
conclusion that the extraction of entanglement from the 
vacuum is a weak and fragile yet intriguing phenomenon. 

A direct experimental verification would first require 
us to consider more realistic models. For example, as a 
reasonable approximation to QED, the detectors could 
be modeled as two-level systems coupled to the electric 
field in the dipole approximation [18j. Another perhaps 
more promising possibility is to use a quantum field ana- 
log such as a linear ion trap [19|, 120] ■ In this context, 
Dirichlet boundary conditions are already effectively im- 
plemented due to the finite number of ions. In addition, 
one could implement a classical potential by introducing 
an external electric field and a thermal state may be ef- 
fectively simulated by immersing the ions in a thermal 
bath. 

It would be interesting to investigate other types of 
boundary conditions such as periodic boundary condi- 
tions. Indeed, this type of boundary conditions could 



easily be simulated with a circular arrangement of ions. 
Furthermore, it may be interesting to explore the effect of 
a classical potential beyond the perturbative treatment. 
This treatment could greatly increase our understanding 
of the modification of the entanglement dynamics caused 
by the potential and allow us to study a greater class of 
potentials. Finally, it should also be interesting to study 
other excited states of the quantum field. For example, 
one could investigate if a 1-particle state ap\0) creates 
more entanglement in the qubits than the vacuum |0) 
and analyze how the negativity depends on the direction 



and magnitude of p. This analysis could be also extended 
to more general multi-particle states. 
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